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Abstract 

The sign structure of correlations of conserved charges are investigated in a QCD like model: the (2 + 1) flavor Polyakov 
Quark Meson model. We compute all susceptibilities of the conserved charges on the (fiB —T) plane up to fourth order and a 
few at higher order as well. By varying the mass of the sigma meson, we are able to study and compare scenarios with as well 
as without a critical point. In the hadron-quark transition regime we identify certain correlations that turn negative unlike 
expectation from ideal hadron resonance gas calculations. The striking feature being that these remain negative deep into the 
hadronic side and thus could be measured in experiments. Measurement of such quantities in heavy ion collision experiments 
can elucidate the location of the QCD transition curve and possibly the critical point. 
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The phase diagram of Quantum Chromodynamics (QCD), the theory of strongly interacting matter, has been a 
subject of intense study both theoretically and experimentally for some time now. The thermodynamic state of 
the strongly interacting medium which is expected to be created in a heavy ion collision (HIC) experiment can be 
specified by four quantities: temperature T and three chemical potentials corresponding to the conserved charges 
baryon number electric charge Q and strangeness S', namely I^q and jis respectively. The QCD degrees of 
freedom are sensitive to these thermodynamic quantities resulting in a rich phase diagram niia. At zero chemical 
potentials = 0), first principle Lattice QCD (LQCD) Monte-Carlo simulations have shown that QCD 

undergoes a smooth analytic crossover transition from the low T hadron resonance gas (HRG) phase to the high T 
phase of quarks and gluons (QGP) O [4]. At non zero but small recent LQGD results O |6] show similar 

behavior. 

However, the QGD phase diagram for iib/T ^ 1 is far from established. Direct first principle techniques of LQGD 
using Monte-Garlo methods fail due to the sign problem mu]. This is where QGD like models which have been tuned 
to reproduce LQGD results at zero could provide valuable insight about the nature of the QGD medium [MH]. 
Model computations at large i^b/T predict the possibility of a first order phase transition [9UT8]. Thus the location 
of the critical point (CP) which is the end point of this first order transition line where there is a second order phase 
transition is an important landmark on the QCD phase diagram. Mapping the phase transition line and the CP is a 
major goal of the heavy ion collision experiments m- 

The CP dynamics gives rise to diverging correlation length that result in non-monotonic variations of some quantities 
which have been proposed as plausible observables to identify the CP [20]. For example, moments of conserved charges 
that can be extracted experimentally through event by event analyses are good candidates to hunt for the CP [2TJ [22] . 
At zero /r^, with three flavors of quarks there are numerous computations of the susceptibilities both on the lattice [23]- 
[26] as well as in models [27H33]. Recently, for non-zero but small susceptibilities have been computed on the 

lattice [341436] . Some of these cumulants have also been computed on the jj^b —T plane in models [29l [37] . 

Such observables are good markers of the CP, as long as they are measured close to the CP, the location of which is 
unknown. Further, in HIC experiments, the produced fireball has finite size and lifetime which can tame the divergence 
of the correlation length and render it finite. This will blur the effects of singularity in the critical region and hence 
diminish the chances of a direct experimental confirmation of the CP [381440] . Thus, rather than looking at the absolute 
values of the susceptibilities, sign structures of the same might be better suited for such studies [40]. Third moments of 
conserved charges like H, Q and energy have been already studied in this regard [40]. These were found to change sign 
at the hadron-quark phase boundary corresponding to peak like structures of second order susceptibilities. Studies 
based on the Polyakov Quark Meson (PQM) model show that higher order cumulants of B and Q become negative 
valued in the transition regime miiig. It has been suggested that such distinctive sign structures follow from the 
scaling functions of the 3-D 0(4) universality class [42]. Higher order generalized susceptibilities for B and their sign 
structure in the phase diagram was studied in the (2 -1- 1) PQM model [37]. For all the above cases, the negative 
regions were found very close to the phase boundary and mostly on the QGP side. It is interesting to note that LQGD 
computations along the chemical freeze-out curve as determined from HRG analysis of yields show that the kurtosis 
of B exhibits a change of sign around a/S'tvat ^ 20 GeV [43]. This has been attributed to proximity to the GP [43] . 
Such sign structures of susceptibilities due to the quark-hadron transition regime and the GP can be observed if the 
chemical freeze-out (GFO) curve also passes very close to the phase boundary or the dynamics is such that the sign 
structures are retained during expansion between the phase boundary and the GFO curve. However, observation of 
negative baryonic kurtosis has remained elusive so far in the Beam Energy Scan (BES) program at RHIG [44]. This 
motivates us to investigate the sign structures of the off-diagonal components of the correlations of conserved charges 
which has so far been ignored. In this paper we work with the (2 + 1) flavor PQM model at the mean field level and 
analyze the sign structures of the cumulants of conserved charges on the {jib ~ T) plane with and without GP. We 
do not find any unique sign structure that could be attributed to the presence of the GP alone. On the other hand, 
there are a few candidates, as summarized in Table [H] that are sensitive to the crossover/transition region whether 
or not there is a GP and show a change of sign. Thus these observables are good indicators of the transition regime. 
Some of these off-diagonal susceptibilities were found to exhibit negative regions that extend deep into the hadronic 
side and hence could be more easily accessible to experiments. 

The rest of this paper is organized as follows: In Section [T| we provide the details of the PQM model and its 
parameters used in this study. We further define the susceptibilities that are studied here and their connections 
with the corresponding moments of the conserved charges that are experimentally measurable. In Section [IT] we first 
compare our model computations with those of lattice at non-zero but small We compare, with lattice, some 

cumulants as well as the values of the strangeness chemical potential obtained by imposing the strangeness neutrality 
condition. We then present our results on various cumulants on the jiB —T plane and comment on their usefulness in 
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mapping out the QCD phase diagram, namely identifying the transition regime and the location of the CP. Finally, 
in Section IIIIl we summarize and conclude. 


I. FORMALISM 

The relevant thermodynamic potential ft (T, /ibi l^s) in the (2 + 1) flavor PQM with the inclusion of the vacuum 
term at a temperature T and chemical potentials l^Q and fis in the mean field approximation is given in Refs. [45l 
1^ . The pressure Pr is given by 


Pr (T, /iB, /iQ, (F, /iB, /iQ, Ms) • 


( 1 ) 


The cumulants of the conserved charges are computed by taking appropriate derivatives of Pr 


{T, 


di+i+>‘{Pr/B) 


( 2 ) 


These generalized susceptibilities are related to the moments of the distribution of the conserved charges such as the 
mean M, variance cr^, skewness S etc. The derivatives in eq. ([^ have been computed numerically using the package 
ADOL-C [Ul US] which allows efficient computation of higher order susceptibilities without further truncation errors. 

When relating to heavy ion collision experiments, there are two constraints to be met. Since the number of partici¬ 
pating nucleons is not fixed a priori, the net baryon number Nb or electric charge Nq can not be fixed independently. 
But, the ratio of Nb to Nq can be fixed to that of the initial heavy ion used in the experiment 2.5). Since the 
incoming heavy ions carry zero net strangeness the condition of vanishing Ns also has to be imposed. 

Nb/Nq = 2.5 (3) 

Ns = 0 (4) 

Eqs. ® and Q fix fig and MS' respectively. In this work we will mainly present our results for 3 different variants 
of the (2 + 1) flavor PQM model: msigdOO, msig400-phys and msigflOO. These differ from each other on the choice of 
the mass of the sigma meson rria- and also the way we treat Ms and mq- We have used rrifj = 400 MeV for msig400 
and msig400-phys while for msigOOO we put rricr = 600 MeV. While msig400 and msig400-phys possess a CP on the 
{fiB —T) plane, msigGOO has no CP on the phase diagram [45]. Thus we are able to study and compare the effects of 
the CP on the sign structures. In order to understand the effects of non-zero Ms and mq, foi* msig400 and msigGOO we 
use Ms = 0 and mq = 0 while for msig400-phys we use Ms and mq as obtained from eqs. ^ and Q. In Table 0 we 
have listed the above scenarios with their descriptions for easy reference 



msig400 

msig400-phys 

msigGOO 

rria (MeV) 

400 

400 

600 


0 

from 

b and (2 


0 


0 

from 

!) and (2 

1) 

0 

CP 

yes 

yes 

no 


TABLE I. The distinct features of the different variants of the PQM model investigated here. 


II. RESULTS 

In our earlier works gSl SS], we had shown good qualitative agreements between LQCD data and PQM model 
predictions at zero chemical potential for a large number of thermodynamic quantities like pressure, entropy density, 
energy density, specific heat, speed of sound etc., and also for some susceptibilities of conserved charges up to sixth 
order. LQCD data are now available also at non-zero but small Mb/F and also for non-zero MS' and mq satisfying 
Eqs. (§ and 0) mm- Before we present our model results for the susceptibilities on the entire (Ib — T plane, 
we perform a comparative study between LQCD and our model at small Mb/F. Here we will particularly look at the 
extracted values of Ms and few ratios of susceptibilities at non-zero Mb/F. 
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FIG. 1. Ratio of fis/I ^b with T/Tc for different values of /is obtained in model and compared to leading order results obtained 
in LQCD [49]. The thick black line denotes the SB limit. 


A. /is//is I Model vs Lattice 

For a realistic estimate of thermodynamic observables related to the QGP experiments, it is important that the 
computations are done for non-zero /i^ and /ig obtained from the conditions as imposed by eqs. ^ and Q. We have 
implemented this and extracted the values of /is and /ig for different T and /i^. In Fig. we show /is normalized to 
liB foi* several values of /i^ = {3, 300,450} MeV. We have compared these model values with the continuum estimates 
obtained in LQCD [49] and find good agreement. jJis/l^B seems to have a monotonically increasing behavior with 
TjTc^ with a faster rise around and finally saturates in the QGP phase to ^ 1/3 (shown by the thick black line in 
Fig. §). This limiting value can be understood easily in the context of a Stefan-Boltzmann (SB) gas of ideal quarks 
and gluons. In the SB limit, only the strange quarks which carry both S and B decide the /is/i^B value. Now in 
order to ensure that Eq. © is obeyed, the fugacity factors for strange and anti-strange quarks should be unity. Since 
strange quarks carry 5 = 1/3 and S = —1 while anti-strange quarks carry B = —1/3 and 5 = 1, it turns out that in 
order to ensure B/ib + Sjis = 0, j^s/i^B must be —B/S = 1/3. We have ignored /ig in this discussion since its value 
is much less than both /i^ and /i^. Thus the SB limit is independent of /i^ and we see this also in the Fig. where 
different /i^ curves all saturate to 1/3 in the QGP side. On the contrary, in the low T regime where the degrees of 
freedom are hadronic, jis!I ^b is much more sensitive to /i^- This is because in the hadronic regime unlike the ideal 
quark gluon gas, strangeness carriers can be both: baryonic (eg.A) and non-baryonic or mesonic (eg. kaons). If we 
had a mesonic gas, /i^ = 0 will always be the solution of Eq. Non-zero values of fis arise only because of the 
strange baryons. Since there is a large mass difference of the order of 700 MeV between the lightest strange meson 
and baryon, depending on T the relative contribution of the strange baryons differ from that of the strange mesons, 
rising monotonically with T. This results in the monotonically increasing behavior of jis/i^B with T/Tc for constant 
fiB- Again, with increasing value of /i^ the strange baryon contribution increases and this results in a larger value of 
l^s/l^B for same T. 


B* Model vs Lattice 

X? xf 


Recently, a lot of effort has been invested in computing susceptibilities at non-zero /i^ in LQCD in order to confront 
them with experimentally measured moments of 5, Q and S [34H36] . Such a comparison will enable a determination 
of the CEO T and fiB, bridging the gap between LQCD and experiments. This is a complementary program to the 
already quite successful endeavor of determining the freeze-out conditions by comparing the hadron yields between 
experiments and hadron resonance gas models [5QH53] . In Refs. [35l|36], the (/i^/E) variation of the ratio 5^ = xf 1X2 
for V = 5, Q have been measured on the lattice. A comparison between theory and experiment of these ratios can 
provide an estimate of /i^ at CEO [35]. In Eig.j^we have plotted these quantities as obtained in PQM and compared 
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FIG. 2. Ratios of susceptibilities: I'xS compared with HotQCD and Rf 2 — Xi 1X2 compared with WB lattice 

data. 


them with the HotQCD lattice data [35] as well as WB lattice data [36|. Both lattice as well as PQM show an almost 
linear variation of R ^2 with (Ib/T in the range 0 < jis/T < 1. The model over predicts the lattice data in the 
entire range. For example, at iib/T ^ 1, the model is about 20% more than the lattice. We should note that there 
is an uncertainty in the determination of Tc adjusting which it is possible to get a better agreement between model 
and lattice. However, here our main intention in performing this comparison is to demonstrate the good qualitative 
agreement between LQCD and PQM model predictions. This gives us faith to trust the PQM results on the /i^ — T 
plane where there is yet to be any lattice data. Thus, having set the platform for a discussion of the PQM results at 
non-zero we go into the next section where we report on the novel qualitative features of the sign structures of the 
generalized susceptibilities whose measurements in heavy ion collision experiments could shed light on the features of 
the QCD phase diagram. 


C. Sign Structures of Susceptibilities: Model 

Having found good qualitative agreement between PQM and LQCD at non-zero /i^, we will now analyze the sign 
structures of various fluctuations and correlations of conserved charges on the {fiB —T) plane. We shall focus mainly 
on the transition regime between the HRG and QGP phases. The transition region can be broadly classified into 
three categories: firstly for small {Ib/T there is a smooth crossover and no true phase transition. Secondly, for large 
enough we expect a first order phase transition. Thirdly, in the intermediate range of one thus expects 

a critical region with a second order phase transition at the CP where the line of first order phase transition meets the 
crossover line. The existence and location of this CP is a topic of intense current research. In PQM, the location of 
the CP is highly sensitive to the value of rria- used. For example, with rria- > 600 MeV one finds no CP [45]. In order 
to understand the distinguishing features of the QCD phase diagram with and without a CP and also the effect of 
non-zero /ig and /i^, we have computed the susceptibilities for three different variants of PQM as outlined in Table 
msigdOO, msig400-phys and msigGOO. We will now report our findings on the sign structure of several susceptibilities 
up to sixth order in these different regions of the QCD phase diagram for the three cases. 

Before we proceed, please note that all figures shown henceforth will be plots on the iib —T plane. In fact, since we 
are interested only in the qualitative features of negative correlations, we scale the T axis (x-axis) by Tc (the cross-over 
temperature at = 0) and the axis (y-axis) by jic (the value where the first order phase transition line meets 
the fiB axis). In all plots we show the smooth cross-over curve in blue dashed lines. The CP is denoted by a red dot 
whereas a black (thick) solid line denotes the first order phase transition line. Contours of various susceptibilities are 
shown with boxed numbers indicating the numerical value of the susceptibilities along those contours. Finally, regions 
of negative correlation are indicated by orange (shaded) regions. 


1. Diagonal Susceptibilities up to 4th order 

We have shown the baryonic susceptibilities up to 4th order in Figure. [^ Experimentally only net proton number is 
obtained as the neutrons are never observed. Cumulants of net proton number act as proxy to the susceptibilities of 
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the net baryon number. Up to 2nd order there are no zero contours and therefore no regions of negative correlation. 
The first appearance of a negative region is for Xs • It spans through the crossover as well as first order phase transition 
regions. For all the three cases studied here, the negative region is located just above the hadron-quark transition on 
the QGP side. For the negative region is only in the crossover region and terminates at the CP. Negative values 
for X 4 1 X 2 hs^ve also been measured on lattice at non-zero {Ib |l3]. For the case of msigbOO, even though the CP is 
absent, there is a negative region in the transition regime that stretches all the way up to very small T. Thus, we 
see that negative regions in Xs and yf only imply the proximity of the hadron-quark transition region that may or 
may not include a CP. Thus negative values of Xs and Xa aot necessarily imply the existence of a CP. We have 
not shown plots of the susceptibilities of Q as the sign structure for the electric charge susceptibilities look similar to 
that of the baryonic ones. The contours of S susceptibilities on the {fiB — T) plane are devoid of any interesting sign 
structures. The yf susceptibilities are shown in Fig. We note that in msig400-phys where the physical conditions 
of eqs. © and Q are imposed, the diagonal strange susceptibilities are affected strongly and are distinctly different 
from what one obtains in msigdOO and msigbOO. In msig400-phys the variation of yf along /ib becomes gentler with 
the contours being almost parallel to the axis. 


2. Mixed susceptibilities of order two 

Neither Xi^ aor —Xif have any interesting sign structure. Neither have any negative regions. On the other 
hand, y^^^ has a narrow region where it turns negative for msig400 and msigbOO. Interestingly, this region of negative 
susceptibilities lies in the hadronic side. Fig.|^shows y^^*^ which has negative values for large fiB close to the transition 
region and within the hadronic phase. For msig400 we see that the region is close to the CP. For msig400-phys we 
observe that regions of negative correlation are in fact reduced and very close to the CP. Since negative correlations are 
also seen for msigbOO, we infer that a measurement of negative values of this correlation does not necessarily confirm 
a CP. 


3. Mixed susceptibilities of order three 

The third order mixed susceptibilities show more interesting behavior. In Figim we show X 12 Xm^- I^^^ 

o 

msig400-phys, there is a large region deep into the HRC phase where X 12 negative. Thus, there is a possibility that 
the CFO curve passes through this region and these negative regions could be observed in experiments. Interestingly, 
for both msig400 and msigbOO where fiQ and jxs are zero this negative region is highly localized on the first order 
line. Xi‘f which is not shown here has small negative valued region along the crossover line from the CP in msig400. 
However, this effect is absent in msig400-phys and msigbOO and therefore probably cannot be used to probe the CP 
in experiment. For Xiu ? variants show negative region in the hadronic side, though they do not extend as 

deep as in the case of Xi 2 ‘ Interestingly, the region of negative correlation extends deeper into the hadronic region 
for msig400-phys For the remaining off diagonal susceptibilities of order 3 we find that none of them have regions of 
negative correlation in the hadronic phase, although some of them have negative regions in the QCP phase close to 
the transition curve. 


4- Mixed susceptibilities of higher order 

In Fig. we show only those fourth order off-diagonal susceptibilities that have regions of negative correlation that 
penetrate the hadronic phase significantly. Xi^^ X 13 have negative values only in the hadronic phase 0 We 
observe that while Xi-^^ extends deep into the hadronic phase, Xis is highly localized near the transition region and 
the CP. Also while the negative region for Xi^ appears extended in msig400-phys as compared to the other cases, 
the opposite effect is observed in Xis where the negative region diminishes slightly. For X 2 n.^ X 31 regions of 

negative correlation exist on either side of the QCD transition line at large /ib and close to the CP. These regions 


^ Other negative regions seen in the plot are deep in the QGP phase and we disregard them. 
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FIG. 3. Contours (black dotted) for diagonal susceptibilities up to fourth order, shown on the ixb —T plane. Also shown are the 
cross-over curve (blue dashed) the CP (red dot) and the first order transition curve (black solid). Boxed number indicate the 
value of the susceptibilities along that contour. Negative regions are shown in orange (shaded). The three columns are for the 
three different variants of the PQM model in Table Column 1: msigdOO , Column 2 : msig400-phys and Column 3: msigGOO 
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FIG. 4. Contours (black dotted) for xf shown on the fis —T plane. Also shown are the cross-over curve (blue dashed) the CP 
(red dot) and the first order transition curve (black solid). Boxed number indicate the value of the susceptibilities along that 
contour. Negative regions are shown in orange (shaded). The three columns are for the three different variants of the PQM 
model in Table Column 1: msigdOO , Column 2 : msig400-phys and Column 3: msigGOO 
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FIG. 5. Contours (black dotted) for Xii ^ shown on the fiB — T plane. Also shown are the cross-over curve (blue dashed) the 
CP (red dot) and the first order transition curve (black solid). Boxed number indicate the value of the susceptibilities along 
that contour. Negative regions are shown in orange (shaded). The three columns are for the three different variants of the PQM 
model in Table Column 1: msig400 , Column 2 : msig400-phys and Column 3: msigGOO 


appear as two separate lobes that merge across the first order transition curve and separate above the CP where the 
nature of the transition is that of a cross-over. X 211 shows a slightly enhanced region of negative correlation for the 
physical conditions of msig400-phys in comparison to the other cases while the negative regions are reduced for Xsi 
in msig400-phys. In Figwe show the contours for X 2 ‘S^ which has a behavior similar to the diagonal susceptibilities, 
with the region of negative correlation closely following the QCD transition curve and proceed into the QGP side. 
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FIG. 6. Contours (black dotted) for ^^^d —X? 2 ^ shown on the {is — T plane. Also shown are the cross-over curve 

(blue dashed) the CP (red dot) and the first order transition curve (black solid). Boxed number indicate the value of the 
susceptibilities along that contour. Negative regions are shown in orange (shaded). The three columns are for the three different 
variants of the PQM model in Table [I] Column 1: msig400 , Column 2 : msig400-phys and Column 3: msigGOO 


D. Tracking the QCD phase transition 

After a brief description of the essential features of the sign structure of various susceptibilities in the discussion 
above, we would like to indicate how a measurement of the same might be useful in locating the QCD phase transition 
curve as well as the CP. In Fig. pj we show along with the phase transition curve, negative regions of a few select 
susceptibilities, namely x^, Xm ? ~X ^2 ^ X^i - We include the diagonal xf susceptibility in this 

plot to accentuate the fact the off-diagonal susceptibilities have regions of negative correlation significantly deeper in 
the hadronic phase. It is our proposal that a systematic measurement of these off-diagonal susceptibilities can provide 
valuable information about the proximity of the phase transition curve as well as the CP. 

In order to illustrate this point, let us consider a hypothetical CFO curve as indicated by the (gray dashed line). 
The curve has been drawn by scaling the phase transition curve down by a factor = 0.93 and hence our CFO curve 
follows the phase transition curve. We have marked four points (labeled: A, B, C, and D) at various values of 
along this curve that indicate possible freeze-out points at which the sign of the susceptibilities can be determined. 
Point A, at high T and low (akin to the experiments at LHC) would not show negative values of the susceptibilities 
under consideration. At lower energies and larger we approach point B where we might obtain the negative values 
for one of the susceptibilities, namely —X? 2 ^- we proceed to point C along the CFO, we now have three negative 

susceptibilities; —X? 2 ^: Xm^ Xm^- Finally at point D, in addition to these three susceptibilities, Xm^ is also 
negative. This indicates that point D is closest to the CP. 
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FIG. 7. Contours (black dotted) for fourth order off-diagonal susceptibilities, shown on the jib —T plane. Also shown are the 
cross-over curve (blue dashed) the CP (red dot) and the first order transition curve (black solid). Boxed number indicate the 
value of the susceptibilities along that contour. Negative regions are shown in orange (shaded). The three columns are for the 
three different variants of the PQM model in Table Column 1: msig400 , Column 2 : msig400-phys and Column 3: msigGOO 
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FIG. 8. Contours (black dotted) ^ shown on the /jlb —T plane. Also shown are the cross-over curve (blue dashed) the CP 

(red dot) and the first order transition curve (black solid). Boxed number indicate the value of the susceptibilities along that 
contour. Negative regions are shown in orange (shaded). The three columns are for the three different variants of the PQM 
model in Table Column 1: msigdOO , Column 2 : msig400-phys and Column 3: msigGOO 



FIG. 9. Contours of negative regions for various susceptibilities for msig400-phys. The gray thick dashed curve corresponds 
to the proxy freeze-out line. Also shown are the cross-over curve (blue dashed) the CP (red dot) and the first order transition 
curve (black solid). 
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- 


/ 

/ 

(0.6,0.8)-E 

0.84 

x^^ 

/ 

/ 

(0.7,0.9)-CP 

0.99 


/ 

/ 

(0.6,0.8)-E 

0.84 

X 22 

/ 

/ 

(0.4,0.9)-CP 

0.99 

X 22 

/ 

X 

(0.6,0.9)-E 

- 

X22 

/ 

X 

(0.9,0.7)-E 

- 

Xzi 

/ 

/ 

(0.2,0.1)-CP 

0.98 


/ 

/ 

(0.4,1.0)-CP 

0.99 


/ 

/ 

(0.8,0.6)-E 

0.96 

x^ 

/ 

/ 

(0.8,0.6)-CP 

0.99 

X31 

/ 

/ 

(0.4,1.0)-CP 

0.99 

Xl 3 

X 

- 

- 

- 

X23 

/ 

X 

(0.4,0.9)-CP 

- 

BQtii 

A .113 

/ 

/ 

(0.8,0.7)-E 

0.92 

XlA 

X 

- 

- 

- 


/ 

/ 

O-cp 

0.96 

X33 

/ 

/ 

(0.2,1.0)-CP 

0.97 

x^T 

/ 

/ 

(0.6,0.8)-E 

0.87 

x^^ 

/ 

/ 

(0.5,0.9)-CP 

0.98 

X24 

/ 

/ 

(0.3,1)-E 

0.99 

A 114 

/ 

/ 

(0.7,0.7)-E 

0.87 

Xl 5 

X 

- 

- 

- 


TABLE IL Summary Table: A list of various susceptibilities with a description of their negative regions. The second column 
indicates the presence of regions of negative correlation with a tick and absence with a cross. The third column indicates 
whether any part of the negative region lies in the hadronic phase. The fourth column corresponds to the approximate range 
of the region (see text for meaning). The last column denotes the extent of the negative region into the hadronic side (see text 
for details). Here R corresponds to a scaling factor as described in the text. 
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One must be wary that the CFO curve might have a more complicated shape than the simple scaled curve that we 
have used. In which case the CFO curve may miss some of the negative regions. Nevertheless, measurement of negative 
values of the susceptibilities we have considered here are indicative of proximity to the phase transition region. Most 
importantly, since these susceptibilities are negative deep in the hadronic phase, the CFO curve is more likely to pass 
through them. 

Before we conclude, we summarize our findings in Table [Ilj In this table we list various susceptibilities and describe 
some of their properties. In the second column we of the table we indicate the presence of region of negative correlation 
by a tick mark and the absence of one by a cross. In the third column we indicate whether there there is a region of 
negative correlation in the hadronic phase. In the last two columns we map out regions of negative correlation. This 
we do by first observing that a common characteristic of negative regions is that they originate at some point of the 
fiB —T plane and follow the phase transition curve. The negative regions terminate either at the CP or follow the 
phase transition curve all the way up to the end (E) at zero T. The fourth column thus indicates the point at which 
the negative region begins {fiB/l^c^T/Tc) and whether it follows the phase transition curve up to the CP or to the 
end (E). Einally, for those susceptibilities that have negative regions in the hadronic phase we give the least value of 
the scale factor (R) for which our hypothetical CEO curve (described above) just touches the negative region. As can 
be gleaned from the table, there are several useful candidates that could be used to locate the QCD phase transition 
curve and CP. 


III. CONCLUSION 


The sign structure of diagonal susceptibilities have been studied previously with negative regions found in the HRG- 
QGP transition regime and mostly in the QGP side. Thus only if CEO occurs close to the phase transition regime are 
these negative regions likely to be observed. HRG analysis of yields reveal that the CEO curve is close to the QCD 
transition curve for small /i^ as in the case of LHC [5QH53] . At SPS and future experiments like EAIR where /i^ is 
large, the CEO curve could be further away from the transition curve and therefore it is likely that these experiments 
may miss the negative regions. In this study we find that for particular choices of off-diagonal susceptibilities as listed 
in Table [llj negative regions could extend deep into the hadronic phase and should be accessible to the experiments 
even at large values of Our study suggest that the correlations of conserved charges exhibit a rich sign structure 
in the fiB — T plane that can be accessed by experiments. These can guide us in our search for the CP as well as 
the hadron-quark phase transition line. Their detection will also signal unambiguously the creation of a novel phase 
of QCD matter unlike the HRG phase. However, we should caution the reader that unlike critical exponents, sign 
structures of cumulants and the region over which they are negative are model dependent. Thus, our results should 
be viewed with caution when we want to extrapolate to real QGD. Another caveat is that , currently this study is 
at the mean field level. As an immediate step, the fate of these negative structures on including beyond mean field 
physics needs to be worked out. We also hope that our work motivates a similar study on the LQGD front on the sign 
structures of the off-diagonal susceptibilities which have so far been ignored. 
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